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High  Resolution  Infrared  Spectroscopy  of 
Carbon  Dioxide  and  Nitrous  Oxide 
at  Elevated  Temperatures 


1.  INTRODUCTION 

■  r  //  b 


C  *  >  r  //  -  IV  r 

J  Although  both/b**v>  and  N^^are  minor  constituents  of  the  terrestrial  atmos¬ 
phere,  they  play  leading  roles  in  several  current  atmospheric  problems.  Both 
CO  and  N  O  are  greenhouse  gases,  and  while  it  is  known  that  the  atmospheric 
concentrations  of  both  gases  is  increasing^-1’  ^the  impact  of  these  trends  on  global 
temperature  is  not  yet  adequately  understood.  3  In  addition  to  being  a  greenhouse 
gas,  NzO  plays  an  indirect  part  in  ozone  chemistry.  ’  _L_  ta¬ 

in  addition  to  the  direct  knowledge  gained  by  studying  the  high  temperature 
spectra  of  gases,  these  studies  help  increase  the  general  understanding  of  the 


(Received  for  publication  2  March  1988) 
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•  ■  Ksies  of  molecules/'  Heating  a  molecule  makes  possible  the  observations  of 
molecular  transitions  originating  from  higiily  excited  rotation- vibration  states, 
these  observations  can  then  be  used  to  characterize  the  shape  of  the  molecular- 
potential  function  at  increasing  distances  from  the  minimum  of  the  potential  function. 

The  infrared  spectra  of  linear  molecules  like  CO.-,  and  NO  are  composed  of 
vibration  bands  made  up  of  a  number  of  nearly  equally  spaced  rotation  lines. ^  At 
high  temperatures,  there  is  a  great  deal  of  overlapping  of  these  bands.  This  over¬ 
lapping  causes  two  main  problems.  First,  the  line  density  in  the  experimental 
spectrum  becomes  very  high,  requiring  high  spectral  resolution.  Second,  this 
overlapping  masks  the  regular  structure  of  each  band,  complicating  line  assign¬ 
ment.  To  assign  lines,  it  is  very  helpful  to  have  wide  continuous  spectral  coverage 
so  that  entire  band  systems  can  be  observed.  A  Fourier  Spectrometer  is  ideally 
suited  to  the  study  of  high  temperature  gases  since  it  meets  both  of  these  needs  by 
providing  high  resolution  over  a  wide  spectral  region.  "•  ■>  • 

The  AFGJ.  high  resolution  Fourier  Spectrometer  has  been  used  to  carrv  out  an 
ongoing  study  of  atmospheric  gases  at  elevated  temperatures.  This  report  sum¬ 
marizes  the  results  that  have  been  obtained  during  the  present  reporting  period, 
spectra  were  recorded  of  X.,0  and  three  different  isotopic  species  of  CO,,  in  several 
different  wavelength  regions.  The  observed  lines  in  the  experimental  spectra  were 
then  identified  with  individual  molecular  transitions  for  as  many  molecular  transi¬ 
tions  as  possible  (over  12,000  transitions).  The  final  step  was  to  use  a  least-square 
fit  to  calculate  new  effective  molecular  constants. 

2.  THEORETICAL  BACKGROUND 

From  a  fundamental  point  of  view,  CO£  and  NgO  possess  similar  physical 
properties.  Both  CO£  and  >^0  are  linear,  triatomic  molecules.  They  are  iso- 
electronic  and  have  a  nuclei  of  nearly  equal  charge  and  mass.  Only  the  electronic 
ground  state  of  either  molecule  is  significantly  populated  in  the  present  experiment, 
yet  its  configuration  plays  an  important  role  in  determining  the  rotation-vibration 
energy  levels.  This  occurs  through  the  coupling  of  the  electronic  spinor  orbital 
angular  momentum  to  the  motion  of  the  nuclei.  However,  in  the  ground  electronic 
state  of  both  CO2  and  NgO  the  net  orbital  and  spin  angular  momenta  of  the  electrons 
are  zero,  so  there  is  no  net  electronic  angular  momentum  to  couple  with  the  motion 
of  the  nuclei.  There  is  a  Fermi  resonance  between  the  bending  mode,  2 ^ »  and 

6.  Bowens-Jenkins,  P.  E.  ,  Cooper,  D.  L.  ,  and  Richards,  W.  G.  (1985)  Ab  initio 
computation  of  molecular  similarity,  J.  Phys.  Chem.  ,  89(No.  11), 


the  symmetric  stretching  mode,  i>  ,  in  both  molecules.  Coriolis  interactions 

1  8 
and  1-type  doubling  also  occur  for  both  molecules.  The  Coriolis  interaction 

occurs  because  as  the  molecule  rotates,  the  asymmetric  stretching  mode  be¬ 
comes  coupled  with  the  bending  mode  v For  vibrational  bands  where  the  1-type 
doubling  interaction  occurs,  each  rotation-vibration  energy  level  splits  into  two 
levels;  an  "e"  level  with  a  symmetric  wave  function,  and  an  "f"  level  with  an  anti¬ 
symmetric  wave  function.  The  degeneracy  of  the  e  and  the  f  energy  levels  is 
removed  by  the  rotation  of  the  molecule. 


In  several  respects,  however,  C02  differs  fi-om  The  main  difference 

between  the  two  molecules  is  the  symmetry.  has  a  symmetric  structure  O-C-O 

where  the  N’  O  structure  is  asvmmetric  N-N-O.  These  symmetry  differences 
greatly  affect  the  nature  of  the  spectra.  Due  to  the  symmetry  of  the  CO^  molecule, 
transitions  involving  the  symmetry  stretch  mode,  u  ^ ,  are  not  dipole  allowed.  In 
addition,  alternating  lines  are  missing  from  rotation-vibration  transitions  (they 
have  zero  statistical  weight).  If  both  of  the  oxygen  atoms  in  the  COg  molecule  are 
not  of  equal  mass,  the  symmetry  of  the  molecule  is  broken  and  the  character  of 


the  spectrum  becomes  more  like  that  of  N„0.  This  study  covered  both  the  sym- 

.  •  .  ■  12_16„  13~16~^  a  12_18„  .  .. 

metric  isotopic  species  C  02>  C  C^.  and  C  O 2  and  the  asymmetric 


1  2  16  „  18  -  13  A6..1S_  ,  13„  16  _  1 7  _ 

species  C  O  O,  t  O  O,  and  COO. 

The  rotation- vibration  energy  term  values,  T(v,  J),  of  a  linear  molecule  can 


be  expressed  as  a  power  series  in  J(.J+  1),  that  is 


T(v,  J)  = 


B  J(J+  1) 

v 


.J(J+  l)]2  +  Hy[J(J+  l)]3  + 


Ly[J(J- 


where  G  ,  B  ,  and  so  on,  are  effective  molecular  constants.  Each  line  in  the 

V  V 

experimental  spectra  corresponds  to  the  transition  between  a  pair  of  rotation- 
vibration  states.  For  those  bands  where  2 -type  doubling  occurs  (2  >  0)  two  sets 
of  effective  molecular  constants  are  used,  one  for  the  e  levels  and  the  other  for 
the  f  levels. 

The  notation  of  the  vibrational  states  that  were  used  for  N„0  and  CO  is  differ- 

z  Q  z 

ent.  For  N^O  the  notation  was  V^V2^V3’  For  t*ie  ^FGL  notation  was  used. 
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3.  He rz berg,  G.  (1945)  Molecular  Spectra  and  Molecular  Structure,  Vol.  II, 
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Absorption  Line  Parameters  Compilation,  AFCRL-TR-73-0096, 


AD  762904. 


In  the  AFGL  notation  the  vibrational  states  arc  identified  by  v1v2'^v3r’  where  "r" 
is  the  ranking  index  assigned  to  each  member  of  a  Fermi  resonating  group  of 
levels.  When  a  state  is  not  involved  in  Fermi  resonance,  r  =  1  and  the  AFGL 
notation  is  essentially  the  same  as  the  notation  used  for  N^O.  When  Fermi  resonance 
is  present  the  ranking  index,  r,  is  appended  to  the  quantum  numbers  of  the  inter¬ 
acting  state  with  the  highest  v..  For  example,  the  AFGL  notation  for  the  two  states 
0  0  1 

10  0  and  02  0,  which  are  highly  mixed  by  Fermi  resonance,  is  10001  and  10002. 

3.  MEASUREMENTS  PERFORMED 


The  Air  Force  Geophysics  Laboratory  high  resolution  interferometer  was 

used  in  conjunction  with  a  high  temperature  absorption  cell  to  make  the  spectral 

measurements.  The  NgO  and  CC^  samples  that  were  used  in  the  study  were 

heated  to  temperatures  up  to  800  K.  The  N90  spectra  were  taken  in  the  8  /am 

region  using  a  sample  of  natural  isotopic  abundance.  The  CO^  spectra  were 

taken  in  the  2.  8  pm  and  4.  3  pm  regions.  Three  different  isotopic  samples  of  CO„ 

13  ^ 

were  used,  a  sample  of  natural  isotopic  abundance,  a  sample  enriched  in  C, 

1 8 

and  one  enriched  in  O.  As  parts  of  this  work  have  been  completed,  the  results 
have  been  incorporated  into  two  previous  AFGL  Technical  Reports.  11  Table  1 
gives  a  summary  of  where  these  results,  including  line  positions,  can  be  found 
for  each  molecule  and  spectral  region.  These  results  have  also  been  published  in 
The  Journal  of  Molecular  Spectroscopy.  ’  in  addition,  the  information  on  line 
position  of  C09  have  been  incorporated  into  the  1986  edition  of  the  AFGL  HIT  RAN 

*  14 

n.  cl“cnlr> "  database. 
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Table  1.  Summary  of  Results  of  High  Temperature  Studies 


Molecule 

Number  of 
Bands 

Wavelength 
Region  (pm) 

Where  Results 
Reported 

14n216o 

18 

8 

This  report 

i2c16o2 

19 

4.3 

AFGL-TR-86-0046 

13c16o2 

1 5 

4 .  3 

A  F  G I  .  -  TR  -  8  6  -  0  04  6 

IW'o 

5 

4.4 

AFGL-TR-86-0046 

12c18o2 

5 

4.3 

AFGL-TR-86-0046 

12c16o18o 

5 

4.3 

AFGL-TR-86-0046 

l2c16ol8o 

5 

4.4 

A  FGL-TR -  87-0016 

^c^o1  ‘o 

i 

4.4 

AFGL-TR-86-0046 

‘V6o2 

1 1 

2.7 

AFGL-TR-86-0046 

l2c16o1Bo 

10 

2.  8 

AFGL-TR-86-0046 

l2c16o2 

2 

2.8 

AFGL-TR-86-0046 

%■ 


* 

* 


* 


* 


* 


ADA  17  3808 
ADA  17  9430 


4.  THE  EXPERIMENTAL  SETUP 

The  experimental  apparatus  consists  of  an  infrared  source,  a  high  tempera¬ 
ture  absorption  cell  and  the  AFGL  high  resolution  interferometer.  A  Nernst  glower 

was  used  as  the  source  of  the  infrared  energy.  The  high  temperature  cell  that  was 

15 

used  has  been  described  previously.  ^  The  central  one-meter  section  of  this  cell 
can  be  heated  to  800 K  and  is  triple  passed  using  the  Pfund  configuration.  The 
total  absorption  path  of  the  cell  is  3.5  meters,  3  meters  of  uniform  high 


15.  Dalton,  W.  S.  ,  and  Sakai,  H.  (1980)  Absorption  cell  for  the  infrared 
spectroscopy  of  heated  gas,  Appl.  Opt.  ,  19:2413. 
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temperat-.  c  and  1/4  meter  Oi.  each  end  of  the  cell  where  the  temperature  drops  to 
nea  r  umoient.  Significant  features  of  the  AFGL  High  Resolution  Inte  rfe  rometer 
Lr.caude  the  use  of  cat's  eye  retro-reflectors,  step  and  integrate  instead  of  con¬ 
tinuous  carriage  motion,  and  a  digital  demodulation  and  integration  scheme. 

The  primary  advantage  of  cat's  eye  retro-reflectors  over  flat  mirrors  is  that 
cat's  eyes  are  insensitive  to  tilt  making  it  much  easier  to  maintain  alignment  as 
the  interferometer  is  scanned.  A  cat's  eye  retro- reflector  also  laterally  displaces 
input  and  output  beams  making  it  possible  to  access  both  output  beams.  The  optical 
signals  from  these  two  beams  are  complementary  and  so  it  is  possible  to  use  two 
detectors  and  operate  them  in  a  push-pull  mode  thus  canceling  out  common  mode 
errors.  Using  dual  detectors  also  helps  to  reduce  the  effects  of  nonlinear 

i  f j 

detectors.  ‘  The  two  Cu:Ge  detectors  of  the  AFGL  high  resolution  inte  rfe  romete  r 
are  mounted  in  the  same  liquid  helium  dewar.  Using  onlv  one  dewar  reduces  cool¬ 
ing  costs  and  increases  convenience.  It  also  helps  match  the  conditions  experienced 
bv  the  two  detectors  making  the  common  mode  rejection  work  better. 

In  our  apparatus  the  infrared  beam  is  chopped  before  entering  the  high  tempera¬ 
ture  absorption  cell.  The  infrared  signals  are  then  detected,  demodulated,  and 
integrated  digitally.  The  digital  data  system  allows  for  fast  settling  time  after  a 
step,  but  long  integration  time  during  data  taking.  It  is  also  used  to  compensate 
for  small  amounts  of  chopper  jitter  and  slight  phase  variations  between  the  two 
complementary  infrared  channels. 

Several  components  of  the  experimental  apparatus  have  been  reworked  during 
this  reporting  period.  10,  11  Previously,  the  maximum  usable  wavelength  of  the 
interferometer  was  about  7  pm,  but  by  installing  a  KBr  beamsplitter  and  Ge:Cu 
detectors  the  usable  wavelength  coverage  has  been  extended  to  approximately 
20  pm.  The  infrared  source  chamber  was  also  totally  rebuilt.  In  addition  to 
these  modifications,  others  are  underway  to  increase  the  accuracy,  reliability, 
and  ease  of  use  of  the  interferometer.  These  additions  include  a  new  stabilized 
reference  laser,  an  improved  KBr  beamsplitter,  a  remotelv  operable  filter  wheel, 
and  a  new  data  system. 


16.  Guelachvili,  G.  (1986)  Distortion  free  inte  rfe  rog  ram  s  in  Fourie  r  transform 
spectroscopy  with  nonlinear  detectors,  Appl.  Opt.  ,  25:4644. 


6 


4.1  Reference  Laser 


A  new  stabilized  reference  laser,  a  Laboratory  for  Science  Model  220,  has 
been  installed  in  the  interferometer.  It  has  a  long-term  frequency  drift  of  less 
than  50  kHz /day.  Long-term  stability  is  particularly  important  for  use  with  a 
step  and  integrate  interferometer  like  the  AFGL  high  resolution  interferometer 
where  each  spectral  scan  can  r®auire  up  to  15  hours.  If  the  potential  stability  of 
the  laser  is  to  be  realized,  the  laser  must  be  maintained  in  a  controlled  environ¬ 
ment.  In  addition  to  controlling  the  physical  environment  for  the  laser,  retro- 
reflections  of  the  laser  beam  must  also  be  controlled. 

The  laser  head  of  the  Laboratory  for  Science  Model  220  laser  is  physically 
separate  from  the  power  supply  and  most  of  the  other  electronics.  It  is  only  the 
laser  head  that  must  be  placed  in  a  controlled  environment.  The  interferometer 
is  operated  in  a  vacuum,  but  the  laser  must  be  kept  at  atmospheric  pressure  so 
the  laser  head  was  placed  in  a  pressurized  enclosure  and  cables  routed  to  the 
exterior  of  the  vacuum  chamber  where  the  rest  of  the  electronics  were  located. 

The  previous  reference  laser  also  needed  to  be  maintained  at  atmospheric  pressure, 
but  due  to  different  mechanical  designs  of  the  two  lasers  it  was  not  possible  to 
place  the  new  laser  in  the  old  laser  enclosure;  thus  it  was  necessary  to  design  and 
construct  a  new  laser  enclosure.  The  design  of  this  enclosure  is  given  in  Figure  1. 

The  long  term  stability  of  the  laser  is  very  sensitive  to  maintaining  the  laser 
head  at  a  fixed  temperature.  To  do  this,  we  maintain  the  laser  head  packaging  at 
an  elevated  temperature.  This  presents  a  problem  because  the  pressure  enclosure 
must  provide  adequate  ventilation  to  prevent  temperature  buildup.  The  laser  head 
can  be  cooled  by  circulating  air  through  the  enclosure  if  extreme  caution  is  used 
to  ensure  that  the  air  flow  is  laminar.  Passing  the  laser  beam  through  a  turbulent 
airflow  would  introduce  fluctuation  in  the  laser  beam  and  turbulent  air  around  the 
laser  head  would  interfere  with  the  operation  of  the  laser  cavity  length  servo. 


(6)  cooling  line. 
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An  additional  complication  arises  since  the  reference  Laser  uses  tr;msverse 
Zeeman  stabilization.  This  method  of  stabilization  is  very  susceptible  to  retro- 
reflected  radiation.  Careful  design  was  needed  to  insure  that  retro- reflected 
energy  from  the  inte rferometer  did  not  reenter  the  laser.  An  interferometer  like 
the  AFGL  High  Resolution  Interferometer  that  uses  cat's  eye  retro- reflectors  is 
very  susceptible  to  backseatte red  radiation  problems.  Backscattered  radiation 
comes  primarily  from  the  cat's  eye  secondarv  where  the  laser  beam  is  brought 
to  a  focus.  A  small  particle  of  dust  or  surface  imperfection  on  the  secondary  can 
easily  scatter  considerable  light  back  into  the  laser.  We  have  found  that  a  neutral 
density  filter  provided  adequate  isolation  from  backscattered  radiation.  Passing 
the  beam  through  a  neutral  density  filter  reduces  backscattered  radiation  since  the 
light  that  is  backscattered  has  to  pass  through  the  filter  twice  while  the  desired 
output  beam  only  passes  through  the  filter  once.  Hence,  a  neutral  density  filter 
with,  an  attenuation  of  10  will  reduce  the  backscattered  radiation  by  a  factor  of  100. 
The  new  laser  is  more  powerful  than  the  old  laser,  so  the  reduction  in  intensity 
by  a  factor  of  10  is  not  a  serious  problem.  If  it  proves  necessary  in  the  future, 
verv  much  higher  levels  of  isolation  can  be  accomplished  by  passing  the  laser 
beam  through  a  polarizer  and  a  quarterwave  plate. 

4.2  Beamsplitter 

The  KBr  beamsplitter  used  in  the  AFGL  High  Resolution  Interferometer  to 
perform  the  N'90  measurements  had  some  deficiencies  and  so  has  now  been 
replaced.  It  had  been  in  storage  for  a  number  of  years  and  had  lost  some  of  its 
flatness,  although  it  performed  satisfactorily  in  the  longer  wavelength  regions 
used  for  the  measurements.  In  addition,  the  beamsplitter  coating  was  such 

that  the  RT  product  of  the  beamsplitter  was  low  in  the  4  (jm  region.  The  coatings 
on  the  old  beamsplitter  consisted  of  a  single  layer  of  germanium.  The  germanium 
coating  was  thicker  on  the  portion  of  the  beamsplitter  used  for  the  infrared  than 
for  the  reference  laser.  The  new  beamsplitter  uses  different  coating  materials 
for  the  two  regions. 

The  most  serious  problem  with  the  old  beamsplitter  was  that  the  germanium 
coating  used  for  the  reference  laser  was  excessively  absorbing.  The  absorption 
was  high  enough  to  make  the  beamsplitter  appear  more  like  a  metallic  than  a 
dielectric  beamsplitter.  The  difference  between  an  interferometer  using  a  metallic 
beamsplitter  and  a  dielectric  beamsplitter  is  the  phase  between  the  two  beams. 

With  a  metallic  beamsplitter  the  outputs  of  the  two  beams  are  in  phase  while  with  a 
dielectric  beamsplitter  they  are  complementary.  1  With  a  dielectric  beamsplitter. 
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since  the  signals  from  the  two  detectors  are  complementary,  it  is  possible  to 
operate  the  electronics  in  a  push-pull  mode  resulting  in  common-mode  rejection. 
With  the  old  KBr  beamsplitter  it  was  possible  to  operate  the  infrared  detectors  in 
a  complementary  manner,  but  not  the  laser  reference  detectors.  With  the  new 
beamsplitter  it  will  be  possible  to  operate  both  the  infrared  and  the  laser  reference 
detectors  in  a  complementary  mode. 

4.3  Filter  Wheel 

Although  a  Michelson  Interferometer  can  cover  a  very  wide  wavelength  range, 
the  photon  noise  produced  at  each  spectral  interval  is  spread  through  the  entire 
spectrum.  Hence,  higher  signal -to -noise  ratios  are  attained  if  the  wavelength 
range  of  the  input  radiation  is  limited  with  an  optical  filter.  Since  this  optical 
filter  is  located  inside  the  interferometer  enclosure,  it  was  necessary  to  deflate 
the  vibration  isolation  pads  and  bring  the  interferometer  enclosure  to  atmospheric 
pressure  to  change  this  optical  filter.  Since  the  AFGL  High  Resolution  Interferom¬ 
eter  is  a  slow  scanning  instrument  (up  to  15  hours),  the  long-term  stability  of  the 
instrument  is  extremely  critical.  More  consistent  results  are  obtained  if  the 
interferometer  is  allowed  to  equilibrate  for  several  days  after  having  been  opened. 
Hence  changing  the  optical  filter  resulted  in  several  days  of  lost  opportunity  to 
take  data.  To  get  around  this  problem  a  six  position  filter  wheel  has  been  installed 
in  the  interferometer  vacuum  enclosure  that  can  be  operated  from  outside  of  the 
vacuum  tank.  The  design  of  this  filter  wheel  is  shown  in  Figure  2. 

4.4  Data  Aquisition  System 

A  new  data  system  is  currently  being  implemented  using  an  IBM  AT  compatible 
computer  to  replace  the  old  system  which  was  based  on  an  outdated  PDP  8/e  com¬ 
puter.  With  the  PDP  8/e  system,  only  data  acquisition  and  control  of  the  inter¬ 
ferometer  were  performed  locally,  and  all  subsequent  processing  of  the  data  was 
performed  using  a  Control  Data  mainframe  computer.  W'ith  the  new  system,  much 
more  of  the  processing  of  the  data  will  be  possible  locally.  The  ability  to  perform 
"quick-look"  checks  on  the  data  before  transferring  it  to  the  mainframe  will  be 
particularly  valuable.  Given  the  rapid  progress  in  the  computer  industry,  in  the 
near  future  it  should  be  possible  to  obtain  a  microcomputer  with  sufficiently  high 
performance  to  perform  the  entire  data  processing  on  a  microcomputer. 


Figure  2.  Design  of  Optical  Filter  Wheel 


The  main  advantages  of  using  an  IBM  PC  for  an  interferometer  controller  are 
low  costs  and  the  flexibility  made  possible  by  the  large  number  of  available  hard¬ 
ware  and  software  options.  However  in  many  ways,  the  PDP  8/e  was  better  suited 
for  use  as  a  programmable  controller  than  an  IBM  PC  computer.  The  PDP  8/e  is 
a  synchronous  machine  with  a  very  simple  non-intrusive  operating  system.  With 
the  PDP  8/e,  the  computer  is  never  servicing  an  interrupt  when  a  time  critical 
operation  is  in  progress.  Also  since  the  PDP  8 /e  is  synchronous,  the  time  re¬ 
quired  to  perform  each  instruction  is  always  the  same,  making  it  possible  to  per¬ 
form  timing  by  using  the  speed  of  the  computer  instructions  themselves.  With  an 
IBM  PC  on  the  other  hand,  since  the  time  required  to  perform  a  given  operation 
is  not  constant,  the  computer  cannot  be  counted  on  to  perform  time  critical  opera¬ 
tions  without  the  aid  of  additional  hardware.  (In  principle  it  is  possible  to  turn  off 


the  interrupt  on  the  IBM  PC,  but  since  the  operating  system  and  the  hardware  of 
the  computer  are  very  closely  integrated  this  is  not  a  very  satisfactory  solution.  ) 
These  disadvantages  are  largely  overcome  by  the  manufacturers  of  the  data 
acquisition  systems  supplying  complete  data  acquisition  subsystems  instead  of 
lust  simple  analog  to  digital  converters.  For  example,  the  MetraByte  Dash- 16 
board  used  in  the  new  data  acquisition  system  includes,  in  addition  to  the  analog 
to  digital  converters,  timing  and  logic  circuits  as  well  as  the  necessary  software 
drives  to  use  them. 

The  new  data  acquisition  and  control  software  is  written  in  FORTRAN  77,  in 
contrast  to  the  old  PDP  8/e  system  which  was  written  in  absolute  loaded  assembly 
language.  Assembly  language  programs  in  general  tend  to  use  the  hardware  more 
efficiently,  but  they  are  also  much  more  dependent  on  the  details  of  the  hardware. 
The  reason  we  have  been  using  a  PDP  8/e  long  after  it  has  become  ooaoi  te  is 
because  of  the  high  manpower  cost  required  to  rewrite  the  data  acquisition  and 
control  program  to  make  use  of  new  hardware.  It  should  be  relatively  easy  to 
modify  the  new  system  to  keep  up  with  advances  in  hardware.  The  ease  of 
programming  in  a  higher  level  language  and  the  availability  of  commercially 
available  graphics  and  mathematical  software,  also  makes  it  practical  to  write 
software  which  is  considerably  easier  to  use  and  much  more  capable. 

5.  N  O  MEASUREMENTS 

2 

Spectra  were  measured  at  several  different  temperatures  and  pressures  to 
facilitate  the  line  assignment  process  and  to  maximize  the  number  of  spectral  lines 
measured  under  optimum  conditions.  After  checking  each  of  these  spectra  for 
consistency,  they  were  co-added  to  obtain  one  spectrum  for  each  temperature  and 
pressure.  The  experimental  conditions  under  which  spectra  were  measured  and 
the  number  of  experimental  spectra  taken  are  listed  in  Table  2.  Measurements  of 
the  spectrum  of  the  empty  absorption  cell  were  interspersed  between  the  N^O 
spectra.  These  empty  cell  measurements  were  used  to  determine  the  100  percent 
transmission  levels.  The  maximum  retardation  of  the  interferometer  for  all 
spectra  was  83  cm,  resulting  in  a  resolution  of  0.  006  cm 
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Additional  spectra  not  included  in  Table  2  were  the  first  obtained  at  800K; 
however,  few  absorption  features  were  observed.  As  these  were  spectra  of  the 
first  samples  run  at  800  K,  it  appears  plausible  that  the  high  temperature  N20  was 
reacting  with  the  walls  of  the  cell.  Residue  from  the  initial  reactions  prevented 
further  loss  of  N^O,  enabling  successful  recording  of  the  interferograms  that 
followed. 

The  spectra  were  calibrated  using  an  internal  calibration  technique.  Toth  has 
1 8 

published  a  paper  in  which  he  reports  the  analysis  of  room  temperature  »20 
spectra  that  were  taken  using  the  high  resolution  Fourier  spectrometer  located  in 
the  Me  Math  solar  telescope  facility  at  the  Kitt  Peak  National  Observatory.  Since 
his  data  were  taken  at  room  temperature  the  lines  that  he  observed  do  not  extend 
to  as  high  rotational  states  (J  values)  as  in  this  work.  The  calibration  was  per¬ 
formed  by  adjusting  the  wavelength  scale  of  the  observed  spectra  until,  on  the 
average,  the  observed  line  positions  for  the  low  J  lines  matched  the  values  obtained 
by  Toth. 

During  the  calibration  process  it  was  noted  that  there  were  some  systematic 
shifts  between  the  positions  of  strong  and  weak  lines  in  spectra  taken  at  high 
temperatures.  This  effect  is  presently  being  investigated,  but  is  probably  due  to 
non-uniform  illumination  of  the  detectors  causing  a  slightly  asymmetric  instrumental 
line  shape.  Evidence  to  support  this  conclusion  is  that  the  quality  of  the  infrared 
beam  has  been  observed  to  be  much  poorer  at  the  higher  temperatures  than  at  room 
temperature.  The  primary  cause  of  this  beam  degradation  was  probably  due  to  the 

18.  Toth,  R.A.  (1986)  Frequencies  of  N„0  in  the  1100  to  1440  cm  1  region, 
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distortion  of  the  mirrors  in  the  absorption  ceil  as  the  cell  was  heated.  Even  for 
the  low  J  lines  of  the  v  fundamental  where  this  effect  was  most  noticeable,  the 
error  was  only  0.0007  cm 

5.1  Treatment  of  the  Data 

The  line  assignments  for  each  band  were  performed  by  starting  at  low  J,  where 
Toth's  molecular  constants  were  valid,  and  working  to  high  J.  After  the  line 
assignments  had  been  made,  data  from  all  the  different  temperatures  and  pressures 
were  combined  into  a  single  data  set  and  a  final  weighted  least-squares -fit  was 
performed  to  obtain  new  effective  molecular  constants.  Although  over  4100  lines 
were  identified  in  the  experimental  spectrum,  only  3454  were  used  in  the  least- 
squares-fit,  due  to  line  merging  problems.  Many  of  the  remaining  lines  were 
slightly  affected  by  the  presence  of  close-bv  spectral  lines.  These  slightly  merged 
lines  were  included  in  the  least-squares  fit,  but  with  reduced  weighting.  Each 
band  was  fit  independently  without  making  any  attempt  to  combine  the  information 
from  the  various  bands  into  a  single  global  self-consistent  set  of  energy  levels  for 
the  molecule. 

In  order  that  this  weighted  least-squares-fitting  procedure  could  be  used,  it 
was  necessary  that  an  estimate  of  the  uncertainty  of  each  experimental  line  be 
made.  The  weight  assigned  each  spectral  line  was  the  reciprocal  of  the  expected 
uncertainty  squared.  The  factors  that  went  into  calculating  the  expected  uncertainty 
of  each  line  were:  the  random  experimental  noise  in  the  spectrum,  line  asymmetry, 
abnormal  width  of  spectral  lines,  and  inconsistencies  of  line  positions  compared 
to  other  lines  in  the  same  band.  The  total  uncertainty  for  each  line  was  defined  as 
the  square  root  of  the  sum  of  squares  of  the  individual  uncertainties.  Further 
details  of  the  methods  used  to  determine  the  weights  for  the  least-squares-fits 
were  described  in  a  previous  report. 

The  most  noticeable  effect  of  using  a  weight  for  each  spectral  line  was  to  sub¬ 
stantially  reduce  the  uncertainty  in  the  spectroscopic  constants  as  predicted  bv  the 
least- squares-fitting  program.  However,  the  spectral  line  positions  calculated 
using  the  resulting  constants  were  found  to  be  quite  insensitive  to  the  values  of 
the  weights  chosen.  This  indicates  that  the  effects  of  line  merging  on  the  position 
of  spectral  lines  were  essentially  random  for  the  high  temperature  spectra  con¬ 
sidered  in  this  study. 


5J2  Data  Analysis  Software 

Even  though  the  data  analysis  software  had  been  used  in  analyzing  asymmetric 
species  of  CC^  that  have  similar  structure  to  N^O  it  was  still  necessary  to  modify 
the  software  to  take  into  account  that  is  less  "harmonic  oscillator  like"  than 
CC>2.  The  energy  expansion  in  terms  of  G,  B,  D,  and  H  doesn't  work  as  well,  so 
it  was  necessary  to  include  an  additional  term,  L.  This  is  partly  due  to  the  fact 
that  Coriolis  perturbation  plays  a  larger  role  with  than  it  does  with  CO^. 

Adding  L's  to  the  fitting  program  made  it  necessary  to  change  the  format  of  the 
molecular  constant  data  files  and  the  formats  of  all  the  subsequent  programs  that 
use  these  molecular  constant  files. 

In  working  with  the  least-squares-fitting  program  it  became  apparent  that 

because  of  numerical  instabilities  adding  L's  was  more  involved  than  just  adding 

an  additional  term  to  the  fitting  function.  The  relative  size  of  the  constants  that 
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were  being  fitted  ranged  from  the  order  of  10  for  G  to  10  for  L.  This  large 
variation  in  the  size  of  the  fitted  parameters  created  problems  with  numerical 
roundoff  in  the  inversion  of  the  matrix  that  was  used  to  obtain  least-squares-fit 
to  molecular  constants.  It  was  found  that  the  numerical  stability  of  the  molecular 
constants  was  considerably  improved  by  scaling  the  constants  and  including  their 
order  of  magnitude  into  the  fitting  equation.  Therefore  in  the  numerical  matrix 
inversion  and  in  the  determination  of  the  constants  all  the  constants  were  of  nearly 
the  same  magnitude.  If  it  would  have  been  necessary,  additional  least-squares 
techniques  could  have  been  used  to  reduce  the  numerical  instabilities  further. 

In  addition  to  changing  the  software  to  include  L's,  several  other  changes  were 
made  to  make  the  software  easier  to  modify  in  the  future.  One  of  these  changes 
was  to  convert  the  programs  from  Control  Data  Corporation  FORTRAN  IV  to 
generic  FORTRAN  77.  Converting  to  FORTRAN  77  also  eliminated  the  inconvenience 
of  being  tied  to  a  single  computer. 

There  are  many  advantages  in  using  more  than  one  computer  to  perform  data 
analysis.  The  computer  work  was  performed  partly  on  an  IBM-PC  compatible  and 
partly  with  a  Control  Data  Corporation  (CDC)  mainframe  computer.  Due  to  the  more 
convenient  operating  environment  and  better  editors,  it  was  found  to  be  more 
productive  to  perform  program  editing  as  well  as  much  of  the  program  development 
on  an  IBM-PC  compatible.  However,  most  of  the  actual  data  analysis  was  per¬ 
formed  using  the  CDC  mainframe  due  to  faster  execution  speeds  and  larger  disk 
storage  size.  An  additional  advantage  of  using  two  computers,  with  different  word 
sizes  and  which  handle  floating  point  calculation  differently,  is  the  ability  to  quickly 
detect  numerical  roundoff  in  the  algorithms  that  are  being  used.  The  same  calcula¬ 
tion  can  be  run  on  both  computers  and  the  results  compared,  thus  checking  for 
numerical  roundoff. 


5.3  Results  and  Discussion 

The  1!>  rotation- vibration  bands  of  X.,0  for  which  molecular  constants  were 
obtained  are  indicated  on  the  energy  level  diagram  of  Figure  1.  The  range  of 
P  and  R  lines  used  in  the  least-squares  fits,  the  total  number  of  'ines,  and  the 
rms  error  for  each  band  is  given  in  Table  3.  The  effective  molecular  constants 
which  were  obtained  are  given  in  Table  4.  The  line  position,  observed  minus 
calculated,  and  expected  uncertainty  of  each  line  used  in  the  least-squares  -fit  are 
given  in  the  appendix.  These  constants  are  effective  molecular  constants  and  so 
should  not  be  expected  to  accurately  represent  the  internal  structure  of  the  \9Q 
molecule.  The  purpose  of  these  effective  constants  is  to  provide  a  means  of 
reproducing,  within  the  experimental  accuracy,  the  position  of  spectral  lines  over 
the  range  of  J  values  covered  by  the  measurements  (see  Ta>  e  3).  There  are 
a  great  many  interactions  between  different  vibrational  states  for  me  .\,,0  molecule. 
The  effects  of  these  inte ractions  are  accounted  for  bv  allowing  the  different 
effective  molecular  constants  to  float  freely  in  the  least-squares-fit  of  each  band. 
Molecular  constants  obtained  in  this  manner  are  not  self-consistent.  For  example, 
the  molecular  constants  obtained  for  the  vibrational  state  1110  from  the 
1  MO  *-  Oi.in  band  a^e  not  consistent  with  those  obtained  from  the  1310  *-  1110  band. 

An  interesting  observation  that  becomes  apparent  from  studying  both  CO,,  and 
N,,C)  is  that  the  interactions  between  levels  not  involved  in  Fermi  resonances  are 
much  stronger  for  NgO  than  for  CO,.  These  interactions  result  in  NgO  being 
less  "harmonic  oscillator  like"  than  €02*  The  effects  of  these  interactions  can  be 
seen  in  Table  4  In  the  magnitude  of  the  inconsistencies  of  effective  molecular  con¬ 
stants  determined  from  different  bands  and  from  the  extremely  large  values  for 
L's  which  are  obtained  for  some  bands. 

Each  band  was  fit  twice,  once  using  L's  and  once  without.  The  spectroscopic 
constants  ( L '  for  the  upper  state  and  L"  for  the  lower  state)  were  included  in  the 
final  least-squares  fit  only  when  their  inclusion  markedly  improved  the  quality  of 
the  fi c  \a.  reduction  in  the  rms  error  of  more  than  20  percent)  and  the  uncertainties 
in  L  were  smaller  than  the  value  of  L  for  both  the  upper  and  the  lower  states. 
Occasionally,  an  exception  was  made  for  bands  where  2 -type  doubling  was  present 
(bands  where  2  >  0).  If  the  e  levels  indicated  the  need  of  an  L  and  the  f  levels 
did  not,  for  consistency  L's  were  used  for  both  sets  of  levels. 
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One  of  the  most  common  wivs  to  deal  '.vith  small  interactions  between  rotation- 

vibration  states  is  through  the  use  of  contact  transformations.  Assuming  that  tins 

tecnnique  is  valid,  the  two  sets  of  spectres*  sue  constants  that  occur  for  bands 
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where  1  -type  doubling  occurs  are  not  independent.  "  For  these  bands,  several 
of  the  spectroscopic  constants  for  the  e  and  the  f  sets  of  levels  should  be  con¬ 
strained  to  be  equal.  This  interdependence  between  sets  of  rotational  constants  is 
a  function  of  1.  When  1  =  1,  the  vibrational  term  values  G  and  should  be 
constrained  to  be  equal.  When!  =  2,  in  addition  to  having  Gfj  equal  to  G,.,  the 

rotational  constants  B  and  B.  should  be  equal.  When  I  =  3,  then  D  =  D,,  and 

e  t  ^  e  i 

so  on.  These  constraints  seemed  to  work  very  well  for  the  CO„  bands  considered 
in  this  study,  but  for  N  ^  ^ J  they  seemed  to  be  causing  some  problems  with  the  fits. 
One  of  th»  things  that  will  be  investigated  as  part  of  future  efforts  is  whether  tnese 
constraints  should  be  dropped  for  N.,u. 

The  residual  to  ttie  least-squares-fits  are  plotted  in  Figures  4-8  for  several 
rotation-vibration  bands.  The  \.,0  line  positions  calculated  using  the  molecular- 
constants  reported  by  other  researchers  are  also  indicated  on  these  same  plots, 
lor  low  lines,  the  N 90  positions  of  ooth  Guelachvili“  and  Toth*°  are  within 
the  experimental  error  of  the  results  presented  in  this  report,  but  at  higher  -I  the 
measurements  start  to  diverge.  This  result  is  not  surprising  since  both  Gueluchvil 
mri  Toth  used  room  temperature  absorption  cells.  t)f  the  two  measurements  Toth's 
is  the  most  recent  and  covers  a  larger  number  of  rotation- vibration  bands.  The 
line  positions  reported  by  Toth  are  also  in  better  agreement  with  the  values 
reported  in  this  study  than  are  Guelachvili's  line  positions.  Although  .'oth's  low 
■  '  values  are  very  accurate,  they  cannot  be  used  to  extrapolate  the  position  of 
the  high  ,J  lines  observed  in  this  studv.  For  example  if  Toth's  constants  were  used 
to  predict  the  position  of  the  high  J  lines,  the  error  would  sometimes  be  as  large 
as  0.05  cm  (see  Figure  8). 


Id.  Amat,  G.  ,  and  N'ielsen,  11.  H.  ( 1 1 )  5  B )  Vibrational  1  -  type  doubling  and  i-tvpe 
resonance  in  linear  polyatomic  molecules,  J.  Mol.  Spectrosc.  ,  2:152. 

20.  Guelaehvili,  G.  (1982)  Absolute  NoO  wavenumbers  between  1118  and  1343  cm 
by  Fourier  transform  spectroscopy,  Gan,  .J.  Phys.  ,  60:1334. 
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Figure  5.  Comparison  for  the  0200  «-  0000  Band 


Figure  6.  Comparison  for  the  1200  *-  0200  Band 
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Figure  8.  Comparison  for  the  1001  *-  0001  Band 


6.  CONCLUSION 


During  the  period  covered  by  this  report  an  extensive  study  has  been  made  on 

C09  and  using  the  AFGL  High  Resolution  Interferometer  in  conjunction  with  a 

high  temperature  absorption  cell.  This  study  includes  the  identification  of  over 

11,  000  C09  lines  belonging  to  78  different  rotation- vibration  bands  and  over 

4100  N90  lines  belonging  to  18  bands.  Many  of  the  high  J  lines  for  these  bands  have 

not  been  observed  previously.  A  weighted  least-squares-fit  technique  was  then 

used  to  obtain  effective  molecular  constants  for  each  of  these  bands.  These 

effective  molecular  constants  predict  the  position  of  spectral  lines  originating 

from  excited  rotational  states  with  an  accuracy  considerably  greater  than  those 

previously  available.  The  CO„  data  have  already  been  incorporated  into  the  1986 
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edition  of  the  AFGL  HITF1AM  molecular  database.  The  N9G  data  as  well  as  the 
C09  data  will  be  incorporated  into  an  AFGL  high  temperature  database  that  is 
presently  being  compiled. 

In  addition  to  providing  data  for  the  AFGL  databases  this  extensive  data  makes 
it  possible  to  come  to  several  other  interesting  conclusions.  Molecular  constants 
obtained  using  room  temperature  CO^  and  \’90  are  not  adequate  for  predicting  the 
position  of  spectral  lines  observed  at  high  temperatures  even  when  the  room 
temperature  measurements  are  extremely  accurate.  The  interaction  between  levels 
not  involved  in  Fermi  resonances  are  much  stronger  for  N90  than  for  C09.  These 
interactions  make  less  "harmonic  oscillator  like"  than  C02. 
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